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ABSTRACT 

For a one-dimensional chain of three nuclear spins (one half) , we make the numerical 
simulation of quantum teleportation of a given state from one end of the chain to 
the other end, taking into account first and second neighbor interactions among the 
spins. It is shown that a well defined teleportation protocol is achieved for a ratio 
of the first to second neighbor interaction coupling constant of J' /J > 0.04. We 
also show that the optimum Rabi's frequency to control the non-resonant effects is 
dominated by the second neighbor interaction coupling parameter (J'). 



1 



1. Introduction 



Quantum teleportation is a technique which allows us to move a quantum state of our 
quantum system from one location to another location using quantum computation 
or quantum information [1]. Quantum computation or quantum information uses 
quantum bits (called qubits) to handle the information on the quantum system. A 
qubit is a superposition of any two quantum basic states of the system (called |0) 
and |1)), * = C o |0) + d\l), such that |C | 2 + | Ci | 2 = 1. One may call |0) and 
|1) the basic qubits of the system. The L-tensorial product of L-basic qubits form 
an L-register of L-qubits, \x) = ■■■ ,io) with ij = 0, 1 for j = 0, . . . , L — 1 

("0" for ground state and "1" for exited state). The set of these L-registers makes 
up the 2^-dimensional Hilbert space where the quantum computer and quantum 
information work. A typical element of this space is ^ = J2 C x \x), where \ C X \ 2 = 

1, and \C X \ 2 gives us the probability of having the state \x) after measurement. 
We are interested in studying the phenomena of quantum teleportation in a solid 
state quantum computer formed by a one-dimensional chain of nuclear spins [2], 
where a preliminary study of this technique was done using two-qubit registers 
and considering first neighbor interaction between spins [3]. In this paper, we will 
consider also the second neighbor interaction among the spins in a chain of three 
nuclear spins system. Thus, the phenomena of quantum teleportation from one end 
of the chain of spins to the other end is studied numerically. Adding a quantum state 
(qubit) to our chain of three nuclear spins means to add another nuclear spin at one 
end, having, then, a quantum computer working with four-qubit registers. We study 
the well performance of the quantum teleportation algorithm through the fidelity 
parameter and determine the minimum value of the second neighbor interaction 
coupling constant to do this. Finally, we also see the modification that the so called 
27r&-method could have with the consideration of second neighbor iteration and its 
implication in our simulation 

2. Equation of Motion 

Consider a one-dimensional chain of four equally spaced nuclear-spins system (spin 
one half) making an angle cos# = l/\/3 with respect the z-component of the mag- 
netic field (chosen in this way to kill the dipole-dipole interaction between spins) 
and having an rf-magnetic field in the transversal plane. The magnetic field is given 
by 

B = (bcos(ujt + <p), — 6sin(u;t + ip), B(z)) , (1) 

where b, u> and <p are the amplitude, the angular frequency and the phase of the 
rf-field, which could be different for different pulses. B{z) is the amplitude of the z- 
component of the magnetic field. Thus, the Hamiltonian of the system up to second 
neighbor interaction is given by 

= - £> k ■ B k - 2Jft £ - 2J'h ]T F k F k+2 , (2) 

fc=0 k=0 k=0 

where // k represents the magnetic moment of the kth-nucleus which is given in terms 
of the nuclear spin as /Lt k = 11^(1?, I k , It ), being 7 the proton gyromagnetic ratio. 
B k represents the magnetic field at the location of the /ct/i-spin. The second term 

2 



at the right side of (2) represents the first neighbor spin interaction, and the third 
term represents the second neighbor spin interaction. J and J' are the coupling 
constants for these interactions. This Hamiltonian can be written in the following 

way 

H = H + W , (3a) 

where Hq and W are given by 

Ho = -h ($> fe /| + 2J(/ 2 If + IIP 2 + /f If ) + 2J'(/ */f + Jf J£ )} (36) 
U=o J 



and 

so JL r . . i 

(3c) 



fc=0 L 



where oj k = r yB{zk) is the Larmore frequency of the kth-spin, 17 = 76 is the Rabi's 
frequency, and 1^ = 1% ± zi| represents the ascend operator (+) or the descend 
operator (-). The Hamiltonian Hq is diagonal on the basis {|*3*2*i*o)}) where ij = 
0, 1 (zero for the ground state and one for the exited state), 

#oK3*2*iio) = E i3i2ilio \i 3 i 2 hio} . (4a) 
The eigenvalues E^i^ are given by 



^=0 fc=o fc=o J 



(46) 



The term (3c) of the Hamiltonian allows to have a single spin transitions on the 
above eigenstates by choosing the proper resonant frequency. 



To solve the Schrodinger equation 



<9* 

fo-fr = # * - (5) 
let us propose a solution of the form 

*(t) = Y / C k (t)\k) , (6) 

k=0 

where we have used decimal notation for the eigenstates in (4a), Ho\k) = E^\k). 
Substituting (6) in (5), multiplying for the bra {m\, and using the orthogonality 
relation (m\k) = S m k, we get the following equation for the coefficients 

15 

ihC m = E m C m + Y,Ck(m\W\k) m = 0,...,15. (7) 

k=0 

Now, using the following transformation 

C m = D m e- iE ^ R , (8) 
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the fast oscillation term E m C m of Eq. (7) is removed (this is equivalent to go- 
ing to the interaction representation), and the following equation is gotten for the 
coefficients D rn 

1 15 

iD m = T E W mk D k e^ , (9a) 
n k=0 

where W m k denotes the matrix elements (m\W\k), and uj m k are defined as 

E m — Ek 

^mk = ^ • (9b) 

Eq. (9a) represents a set of sixteen real coupling ordinary differential equations 
which can be solved numerically, and where W^s are given by 

(W) = -™*{0, V'or" z*} , (9c) 

where z is defined as z = e l ( tJ *+^) ) an d z* is its complex conjugated. 

A pulse of time length r, phase ip, and resonant frequency u = u^j will be denoted 
by Rij(Q,T,ip), and it is understood that such a pulse implies the solution of Eq. 
(9a) with given initial conditions. 



3. Quantum teleportation and numerical simulation 



The basic idea of quantum teleportation [1] is that Alice (left end qubit in our chain 
of three qubits) and Bob (the other end qubit) are related through and entangled 
state, 

(10) 



^ e = ^(J0a00 b ) + |U01 b )) . 



So, neglecting decoherence effects which could destroy this entangled state, we adjoin 
to Alice an arbitrary state ^> x , 



* x = Cff|0) + Cf|l) , 
resulting the quantum state (v&i = \P X <8> * e ) 



V2 



(11) 



(12) 



C£|0000) + Cft|0101> + Cf|1000> + Cf|1101) 

Then, one applies a Controlled-Not (CN) operation in this added stated and Alice, 

CN 2 3\i3,i2,k,io) = \k, h ®h,i2,io) , (13) 
where i 2 © h = h + h{mod 2), resulting the state (^2 = CNzz^i) 
1 



^ 2 = c^|oooo) + c$\om) + cf 1 1100) + cf |iooi) . (u) 

V2 L 

Finally, one applies a superposition (Hadamar) operation in the added state location, 

r |oooo) |oooo) + 1 1000) 

3 \ 1 1000) ~ y/2 I |0000) - 1 1000) ' { b) 
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resulting, after some rearrangements, the state (^3 = -A3 ^2) 

*3 = U|00)®|0)fe|0)+(^|l))+|01)(8)|0)®fe|l)+C7f|0) 



2 

+|io) ® |o) ® ^|o) - cf + |ii) ® |o) ® (q?|i) - cf |o) 

(16) 

When Alice measures both of her qubits, there are four possible cases (|00), |01), 
1 10) , and 1 1 1) ) , and for each case Bob will get the original ^> x state by applying a 
proper operation to his state: identity, Not (JV), a z , or a z N, where the action of 
these operators are as follows: N\Q) = |1), N\l) = |0), cr 2 |0) = |0), and a z \l) = -|1). 
The most important features in this algorithm are that Alice and Bob do not need 
to know the state *f/ x , and there are always four possibilities for Alice measurement. 

Now, the algorithm we need to implement this quantum teleportation technique 
in our one-dimensional three nuclear spins system varies a little from the scheme 
presented above since our three-qubits quantum computer will start from the ground 
state 

^00 = 1 000) . (17a) 

Therefore, one needs first to attach the unknown state (11) to (17a) to get the initial 
four-spins wave function (^0 = <8> ^00)) 

ijr = Q0000) + Cf|1000) , (176) 

where \C§ | 2 + |Cf | 2 = 1. The entangled state (12) is obtained with the following 
three pulses 

*i = i?4, 5 (7r, -37r/2)i? 8jl2 (7r/2, -ir/2)Ro A (ir/2, -ir/2)V . (16c) 
The Controlled-Not operation CN 23 is gotten through the following two pulses 

*2 = 129,13 (7T, 3tt/2)R 8A2 (7T, -37T/2)*! . (17d) 

Finally, the wave function (16) is gotten after the application of the following two 
pulses to the above wave function 

*3 = i?4,i 2 (vr/2, -37r/2)fli i9 (ir/2, -3vr/2)^ 2 • (17e) 

To make the numerical simulation of this algorithm, we have chosen the following 
parameters in units of 27rxMhz, 

uj = 100 , wi = 200 , lo 2 = 400 , uj 3 = 800, J = 10 , J' = 0.4 , O = 0.1 . (18a) 

These parameters were chosen in this way to have a clear definition on our Zeeman 
spectrum and transitions among them. Of course, our main results are applied to 
the actual current design parameters of reference [4]. On the other hand, we selected 
the state (11) with the following coefficients 

CS = l and <?f = X- (186) 
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Fig. 1 shows the Zeeman spectrum of the four nuclear spins system with the tran- 
sitions used during our teleportation simulation. Fig. 2a shows the probabilities 
| Co | 2 , | C5 1 2 , I C$ 1 2 and | C 13 1 2 during the first three pulses, where the formation of 
the entangled wave function (12) is shown, from the initial state (17b). Fig. 2b 
shows the probabilities | Co | 2 , | C5 1 2 , |Cg| 2 and | C12 1 2 during the following two pulses 
to get at the end the wave function (14). Fig. 2c shows the probabilities | Co | 2 , | Cs | 2 , 
I C5 1 2 , I C13 1 2 , I C12 1 2 , I C4 1 2 , I C9 1 2 and I Ci 1 2 during the last two pulses to get at the 
end the desired function (16). Note that at the end of our algorithm one expects 
the following values for these probabilities 

|C | 2 = |C 8 | 2 = |C 5 | 2 = |C 13 | 2 = ^, (19a) 

and 

|C 12 | 2 = |C 4 | 2 = |C 9 | 2 = |d| 2 = l^i! (196) 

To get a better feeling what is going on on each process, we calculate the z- 
component of the expected values of the spin for each qubit. These expected values 
are given by 

<^> = ^E(- 1 )*i c *(*)i 2 . ( 20fl ) 

k=0 

{ID = ^||Co| 2 + |C 1 | 2 -|C 2 | 2 -|C 3 | 2 + |C4| 2 + |C5| 2 -|C 6 | 2 -|C 7 | 2 

+ |C 8 | 2 + |C 9 | 2 - |do| 2 - |di| 2 + |C 12 | 2 + |C 13 | 2 - I C 14 1 2 - |C 15 | 2 1 , 



(206) 



{ID = i||Co| 2 + |d| 2 + |C 2 | 2 + |C 3 | 2 - |C 4 | 2 - |C 5 | 2 - |C 6 | 2 - |C 7 | 2 

+ |C 8 | 2 + |C 9 | 2 + Idol 2 + |di| 2 - |d 2 | 2 - Idsl 2 - |d 4 | 2 - |ds| 



(20c) 



and 



{id = \ii\c k \ 2 -\f:\c k ? . (20d) 

Z k=0 Z fc=8 

Figs. 3a, 3b, and 3c show these expected values during entangled formation, controlled- 
not operation, and final teleportation. As one can see, these behavior is what one 
could expected for each process. Fig. 4a shows the probabilities (19a) and (19b) at 
the end of the teleportation (wave function (16)), and Fig. 4b shows the probabilities 
of the non-resonant states involved in the dynamics. 

To see the values of the second neighbor interaction coupling parameter from which 
one could have a well defined teleportation algorithm, we study the fidelity param- 
eter [5], 

F = <* expected]®) , (21) 
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where ^ expected is the ideal wave function (16). Fig. 5 shows this fidelity parameter 
as a function of the ratio of second to first neighbor coupling parameters (J' / J). As 
one can see, a well defined teleportation algorithm is gotten for J' / J > 0.04. 

4. Quantum teleportation and the 27rfc-method 

One of the important results from the consideration of first neighbor interaction 
and the selection of the parameter as J/Auj <C 1 and 0,/Alu <C 1 is the possibility 
to choose the Rabi's frequency f2 in such a way that the non-resonant effects are 
eliminated. This procedure is called the 27rfc-method [4], and this Rabi's frequency 
is chosen as = |A|/\/4A; 2 — 1 for a -zr-pulse, where k is an integer number, A is the 
detuning parameter (A = (E p — E m )/fi — uj) between the states \p) and \m) when 
the resonant frequency is uj, and this detuning parameter is proportional to the first 
neighbor coupling constant J. Let us see how this detuning parameter could be 
modified due to second neighbor interaction. Assuming that the states \p) and \m) 
are are the only ones involved in the dynamics, from Eq. (9a), one has 

h h 

Thus, given the initial conditions D p (0) = C p (0) and D m (0) = C m (0), the solution 
is readily given by 

and 

where fi e is defined as Cl e = \fVL 2 + A 2 . For a 7r-pulse (t = r = 7r/0), one can 
select the term f2 e 7r/2f2 to be equal to any multiple of ir, O e 7r/2r2 = kir, to get the 
condition 17 = \A\/\/4k 2 — 1. this condition gets rid of the non-resonant terms since 
from Eqs. (20a) and (20b) one gets 

D p (r) = (-l) k C p (0)e iA ^ 2n , and D m {r) = (-l) k C m (0)e~ iA ^ 2n . 

For a 7r/2-pulse (t = t = 7r/20), one can select the term J7 e 7r/4r2 to be equal to 
any multiple of 7r, $l e 7r/4ri = kir, to get the condition Q = |A|/\/l6fc 2 — 1. this 
condition gets rid also of the non-resonant terms. 

Now, if for example one selects a resonant transition which contains the Larmorc 
frequency u>o, these frequencies could be ojq + J + J', ujq — J + J', luq — J — J' or 
ooq + J — J' which correspond to the transitions (decimal notation) |0) <-> |1) (|10) <-> 
|11)), |2) |3), |6) |7), and |2) <-> |3). So, all of these states are pertubated, 
and the frequency difference A may have the values 2 J, 2 J', 2 J + 2 J', or 2 J — 2 J'. 
For other Larmore frequencies the additional values of the detuning parameter are 
4 J, and 4J + 2 J'. Thus, let us denote by the Rabi's frequency selected by this 
method, 

Qa = 7W=i ' n ' puhe (24a) 



D p (t) = C p (0) 



Opt 



cos 



A tt e f 



An(t) = C m (0) 



cos ■ 



a n c f 

z— - sm 

17. 2 
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and 

fiW = 7f J|= T , -/2- P ulse (246) 

where A can have the values 4J+2J', 4J, 2J+2J', 2 J or 2J'. To see the dependence 
of our teleportation algorithm with respect the Rabi's frequency, we use again the 
fidelity parameter parameter (21). With the same values for our parameters as (18a) 
but £2, Fig. 6 shows the fidelity parameter as a function of the Rabi's frequency. 
Dashed vertical lines mark the omega values where the peaks ocurres. These peaks 
correspond to some specific omega defined through the 27rfc-method. For example, 
the line (1), (2) and (3) correspond to the following Rabi's frequency values 

^4J+2J' ~ ^4 J ~ ^4j+2J' ~ ^4 J ^ ~ ^4 J' = 0.10079 , 



^4J+2J' ~ ^4J ' ~ ^2J+2J' ~ ^2J ' ~ ^2J' ~~ 0.066889 , 



and 



,(305) r 
■4J+2J' r 


~ o( 304 ) 


~ 0( 156 ) r 

~ "2J+2J' ' 


~ o( 15 °) 

- S2 2J 


,(407) p 
'4J+2J' r 


~ O ( 40 °) 


~ o( 208 ) r 

~ S2 2J+2J' ' 


^ O (200) 
- Si 2J 



T I O 7/ ~ ^/i 7 ^ ~ 7 I O 7/ ~ 7 ^ ~ ^2J' 0.050098 . 



As one can see in Fig. 7a, where we have plotted (for the detuning values 
mentioned above) and where the corresponding dashed vertical lines of Fig. 6 have 
been drawn, around these lines there are several other values of 0^ +2 j,, fi^j, 

^4j+2J' anc ^ ^2J which, in principle, should cause a peak in the fidelity parameter 
(because they belong to the 27r/c-method) . However, they do not appear at all on Fig. 
6. This means that the peaks values on the fidelity parameter are fully dominated 
by the second neighbor coupling interaction parameter (</')• On the other hand, the 

reason why only even numbers of k appears for these peaks for Q 2 k ], can be seen in 

(k) ~(k) 

Fig. 7b where Vt 2 j, and Q 2 j, have been plotted as a function of k. On this plot one 
sees that the jth-dashed lines correspond to Q 2 j, = ti^j'- Therefore, the peaks on 
fidelity correspond to 27r/c-method dominated by J' and by 7r/2-pulses. This seems 
reasonable since our teleportation algorithm starts with 7r/2-pulses and finishes with 
7r/2-pulses. 
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5. Conclusions 



We have made a numerical simulation of teleportation in a solid state quantum 
computer modeled by one-dimensional chain of three nuclear spins (one half), and 
considering first and second neighbor interactions. We have shown that a good tele- 
portation algorithm can be gotten if the ratio of second to first neighbor interaction 
constants is chosen such that J'/ J > 0.04. We also studied the effect of the sec- 
ond neighbor interaction on the detuning parameter which is used in the so called 
27rfe-method to eliminate non-resonant transitions, and we have shown that the ap- 
plication of this method in our teleportation algorithm is not so simple since the 
detuning parameter varies with both parameters J and J' (first and second neighbor 
coupling interactions). However, the peaks on the fidelity parameter are dominated 
by the second neighbor coupling parameter and the 7r/2-pulses. 
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Figure Captions 

Fig. 1 Energy levels and resonant frequencies used within the algorithm. 

Fig. 2 Probabilities \Ck\ 2 ■ (k). (a) Formation of the entangled state, wave function 
(12). (b) Formation of the wave function (14). (c) Formation the wave function 
(16). 

Fig. 3 Expected values {If,): (k=0,l,2,3). (a) During formation of wave function 
(12). (b) During formation of the wave function (14). (c) During formation of the 
wave function (16). 

Fig. 4 Probabilities |C fc | 2 . (a) For the expected registers k = 0,8,5,13,12,4,9,1. 
(b) For the non-resonant states k = 2, 3, 6, 7, 10, 11, 14, 15. 

Fig. 5 Fidelity parameter MS M function of J' /J. 

Fig. 6 Fidelity parameter function of S7. 

Fig. 7 (a): Rabi frequency as a function of k for A = 4 J + 2 J' [1], A = 4 J 
[2], A = 2 J + 2 J' [3], A = 2 J [4], A = 2 J' [5]. Dashed lines (j) for j = 1, ... ,8 
correspond to Fig. 6. (b): Rabi frequency fi^ [1] and [2] as a function of k. 
Dashed lines (j) for j = 1, . . . , 8 correspond to Fig. 6. 
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